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Abstract 

We introduce a fractional theory of the calculus of variations for multiple integrals. Our 
approach uses the recent notions of Riemann-Liouville fractional derivatives and integrals in 
the sense of Jumarie. Main results provide fractional versions of the theorems of Green and 
Gauss, fractional Euler-Lagrange equations, and fractional natural boundary conditions. As 
an application we discuss the fractional equation of motion of a vibrating string. 
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1 Introduction 

The Fractional Calculus (FC) is one of the most interdisciplinary fields of mathematics, with many 
applications in physics and engineering. The history of FC goes back more than three centuries, 
when in 1695 the derivative of order a = 1/2 was described by Leibniz. Since then, many different 
forms of fractional operators were introduced: the Grunwald-Letnikov, Riemann-Liouville, Riesz, 
and Caputo fractional derivatives [3TJ [35] , and the more recent notions of Klimek 32; , Cresson 
[H] , and Jumarie [H (23 [571 [5S] . 

FC is nowadays the realm of physicists and mathematicians, who investigate the usefulness 
of such non-integer order derivatives and integrals in different areas of physics and mathematics 
[T2l 123] l3T] . It is a successful tool for describing complex quantum field dynamical systems, 
dissipation, and long-range phenomena that cannot be well illustrated using ordinary differential 
and integral operators [501 HS1 13H I3H] • Applications of FC are found, e.g., in classical and quantum 
mechanics, field theories, variational calculus, and optimal control [191 [22l [27] . 

Although FC is an old mathematical discipline, the fractional vector calculus is at the very 
beginning. We mention the recent paper [32], where some fractional versions of the classical results 
of Green, Stokes and Gauss are obtained via Riemann-Liouville and Caputo fractional operators. 

*Accepted for publication in the Journal of Mathematical Physics (14/January/2010). 
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For the purposes of a multidimensional Fractional Calculus of Variations (FCV), the Jumarie 
fractional integral and derivative seems, however, to be more appropriate. 

The FCV started in 1996 with the work of Riewe [55]. Riewe formulated the problem of 
the calculus of variations with fractional derivatives and obtained the respective Euler-Lagrange 
equations, combining both conservative and nonconservative cases. Nowadays the FCV is a sub- 
ject under strong research. Different definitions for fractional derivatives and integrals are used, 
depending on the purpose under study. Investigations cover problems depending on Riemann- 
Liouville fractional derivatives (see, e.g., [20l [21] ) , the Caputo fractional derivative (see, e.g., 
[3 [510), the symmetric fractional derivative (see, e.g., [32 ), the Jumarie fractional derivative 
(see, e.g., [5] [27]), and others [3j HU [19]. For applications of the fractional calculus of variations 
we refer the reader to jTH H3 HOI I2Z1 E21 ESI EZ1 SO]- Although the literature of FCV is already 
vast, much remains to be done. 

Knowing the importance and relevance of multidimensional problems of the calculus of varia- 
tions in physics and engineering [43) . it is at a first view surprising that a multidimensional FCV 
is a completely open research area. We are only aware of some preliminary results presented in 
[201 . where it is claimed that an appropriate fractional variational theory involving multiple in- 
tegrals would have important consequences in mechanical problems involving dissipative systems 
with infinitely many degrees of freedom, but where a formal theory for that is missing. There is, 
however, a good reason for such omission in the literature: most of the best well-known fractional 
operators are not suitable for a generalization of the FCV to the multidimensional case, due to 
lack of good properties, e.g., an appropriate Leibniz rule. 

The main aim of the present work is to introduce a fractional calculus of variations for multiple 
integrals. For that we make use of the recent Jumarie fractional integral and derivative [24, 25 ( 28 , 
extending such notions to the multidimensional case. The main advantage of using Jumarie's 
approach lies in the following facts: the Leibniz rule for the Jumarie fractional derivative is 
equal to the standard one and, as we show, the fractional generalization of some fundamental 
multidimensional theorems of calculus is possible. We mention that Jumarie's approach is also 
useful for the one-dimensional FCV, as recently shown in [5] (see also [26]). 

The plan of the paper is as follows. In Section [2] some basic formulas of Jumarie's fractional 
calculus are briefly reviewed. Then, in Section [5] the differential and integral vector operators are 
introduced, and fractional Green's, Gauss's and Stokes' theorems formulated. Section[4]is devoted 
to the study of problems of fractional calculus of variations with multiple integrals. Our main 
results provide Euler-Lagrange necessary optimality type conditions for such problems (Theorems[3j 
and[5]) as well as natural boundary conditions (Theorem [4]) . We end with Section [5] of applications 
and future perspectives. 



2 Preliminaries 

For an introduction to the classical fractional calculus we refer the reader to [31] 34, 35, 39 . In this 
section we briefly review the main notions and results from the recent fractional calculus proposed 
by Jumarie (25] [27j [28] . Let / : [0, 1] -> R be a continuous function and a € (0, 1). The Jumarie 
fractional derivative of / may be defined by 

f (a H*) = W ^ ) il X (* -t)-(/(t) -/(0))*. (1) 

One can obtain (p} as a consequence of a more basic definition, a local one, in terms of a fractional 
finite difference (cf. equation (2.2) of [28]): 

fia){x) = iTo h (t) ' + (° - k ^ ■ 

k=0 V ' 

Note that the Jumarie and the Riemann-Liouville fractional derivatives are equal if /(0) = 0. The 
advantage of definition (p} with respect to the classical definition of Riemann-Liouville is that the 
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fractional derivative of a constant is now zero, as desired. An anti-derivative of /, called the (dt) a 
integral of /, is defined by 



f{t){dt) a = a / (x-t) a ~ l f(t)dt. 
Jo 

The following equalities can be considered as fractional counterparts of the first and the second 
fundamental theorems of calculus and can be found in [37J [35] : 

£^£f(t)(dtr = a\f(x), (2) 

' f^(t)(dt) a = a\(f(x)-f(0)) ) (3) 

where a! := T(l + a). The Leibniz rule for the Jumarie fractional derivative is equal to the 
standard one: 

(f(x)g(x))^ = f^(x)g(x) + f(x)g^(x). 

Here we see another advantage of derivative (U} : the fractional derivative of a product is not an 
infinite sum, in opposite to the Leibniz rule for the Riemann-Liouville fractional derivative I35| 
P. 91]. 

One can easily generalize the previous definitions and results for functions with a domain [a. b]: 
and 

' X f(t)(dt) a =a { X {x-t) a - l f{t)dt. 



3 Fractional integral theorems 

In this section we introduce some useful fractional integral and fractional differential operators. 
With them we prove fractional versions of the integral theorems of Green, Gauss, and Stokes. 
Throughout the text we assume that all integrals and derivatives exist. 

3.1 Fractional operators 

Let us consider a continuous function / = / [xx, • ■ • , x n ) defined on R — n™ =1 [a^, bi] C K n . Let 
us extend Jumarie's fractional derivative and the (dt) a integral to functions with n variables. For 
Xi € [cii, bi], i = 1, . . . , n, and a s (0, 1), we define the fractional integral operator as 

aJSM =<* I \x l ~t) a - 1 dt. 
J ai 

These operators act on / in the following way: 

ai-CM/^i* •■•,£«) = a / f(xi,.--,Xi-x,t,x i+ i,...,x n )(xi -t) a ~ x dt, i = l,...,n. 
J at 

Let S = {kx, ■ ■ ■ , k s } be an arbitrary nonempty subset of {1, . . . , n}. We define the fractional 
multiple integral operator over the region R~, = IT| =1 [a^ , x^J by 

Ir s felt ■ ■ ■ 1 ^s] = a kl I Xkl fei] ■ ■ ■ a ks I Xks [k s ] 

{x kl - tk^ 1 ■ ■ ■ (x ks - tkj 01 - 1 dt ks . ..dt kl 
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which acts on / by 

Igg [h,... , fc s ]/(a;i, . ..,x n ) 

/(a, • • • ,Zn)(x kl - tkX" 1 ■ ■ ■ {xk s - t ks ) a ^ dt K ■ • • dt kl , 

where £j — tj if j G H, and £j = Xj if j ^ H, j = 1, . . . , n. The fractional volume integral of / over 
the whole domain i? is given by 

iy = oc n r ■■■ [ n f(h,...,t n )(b 1 -t 1 ) a - i ---(b n -t n ) a - i dt n ...dt 1 . 

The fractional partial derivative operator with respect to the ith variable Xi, i = 1, . . . of 
order a € (0, 1) is defined as follows: 

a i D% i [l\= 1 A [ X \x t ~t)- a dt, 

which act on / by 



a % DZ z \l\f{xi,...,X n ) 

1 9 



T(l - a) dxi J ai 



yX-i t) \f y*X\ , . . . 7 Xi — i : t 7 X^i 7 . . . , X n J f {x \ , . . . , 3^2— l , &i , X{j^ \ , . . . , x n )] dt . 



i = 1, . . . , n. 

We observe that the Jumarie fractional integral and the Jumarie fractional derivative can be 
obtained putting n = 1: 

a i«ii}f(x) = a [\x-tr- i f(t)dt= f f(t)(dtr 

J a J a 

and 

a Dz\i]f(x) = Yf^-^ Jj x - *r a (/(*) - /(«)) * = / (Q ° (*)• 

Using these notations, formulae (J2j — d3j) can be presented as 

a D a x [l]J«[l]f{x)=a\f{x) 

a IZ[l] a D«[l}f(x) = a\(f(x)-f(a)). (4) 
In the two dimensional case we define the fractional line integral on dR, R = [a, b] x [c, d], by 

where 

I% R [l]f = alZ[l][f(b,c)-f(b,d)] 

[f(t, c )-f(t,d)} (b-tr-ut 



and 

Ig R {2}f = cimif(b,d)-f(a,d)} 

= a [ [f(b,t)~ /(a,*)]^-*) - 1 *. 
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3.2 Fractional differential vector operations 

Let Wx = [a,x] x [c,y] x [e,z], W — [a,b] x [c,d] x [e, /], and denote {x\,X2,xz) by (x,y,z). We 
introduce the fractional nabla operator by 

VS, X =i 8 D«[l]+j c D°[2] + * e I>?[3] ) 

where the i, j, k define a hxed right-handed orthonormal basis. If / : R 3 — > R is a continuous 
function, then we define its fractional gradient as 

Grad^ x / = V^ x / = i a D°[l}f(x,y,z) + j c D*[2]f(x,y, z) + k e D«{3]f(x,y,z) . 

If F = [F x , F y , F z ] : R 3 — > R 3 is a continuous vector field, then we define its fractional divergence 
and fractional curl by 

Div^ x F = o F = a D«[l]F x (x, y,z) + c D<* y [2}F y (x, y, z) + e D« [3]F Z (x,y,z) 

and 

Curl^ x F = V^ x x F = i ( C D% [2]F Z (x, y, z) - e D a z [3}F y (x, y, z)) 

+ j ( e D«[3]F x (x, y, z) - a D%[l]F z (x, y, z)) + k ( a D%[l]F y (x, y, z) - c D°[2]F x (x, y, z)) . 

Note that these fractional differential operators are non-local. Therefore, the fractional gradient, 
divergence, and curl, depend on the region Wx- 

For F : R 3 — > R 3 and /, g : R 3 — > R it is easy to check the following relations: 

(i) Div^ x (fF) = /Div^ x F + Fo Grad^ x /, 

(ii) Curl^ x (Grad^ x /) = [0,0,0], 
(hi) Div^ x (Curl^ x F) = 0, 

(iv) Grad^ x (fg) = ,gGrad^ x / + /Grad^ x g, 

(v) Div^ x (Grad^ x /) = a D^[l] a D^[l]f + C D«[2] C I>«[2]/ + e D«[3] e D«[3]f- 

Let us recall that in general (£>^ x ) 2 # D$ x (see [25]). 

A fractional flux of the vector field F across dW is a fractional oriented surface integral of the 
field such that 

(Ig w ,F) = I% w [2, 3]F x (x, y, z) + I% w [l, 3]F y (x, y, z) + J£ w [l, 2}F z (x, y, z) , 

where 

I% w [1 , 2}f(x, y,z) = a Ib [l]dj [2] [f(b, d, /) - f(b, d, e)} , 
IS W [1, 3]f(x, y, z) = i?[l] ^[3] [f(b, d, f) - f(b, c, /)] , 

and 

I§ w [2,3]f(x,y,z) = c I%[2] e If[3][f(b,dJ) - f(a,d,f)} . 

3.3 Fractional theorems of Green, Gauss, and Stokes 

We now formulate the fractional formulae of Green, Gauss, and Stokes. Analogous results via 
Caputo fractional derivatives and Riemann-Liouville fractional integrals were obtained by Tarasov 

in gg. 

Theorem 1 (Fractional Green's theorem for a rectangle). Let f and g be two continuous functions 
whose domains contain R = [a, b] x [c, d] C R 2 . Then, 

IsrW + Ig R [2]9 = ^Ir [aX>?[l]9 - cD a d [2]f] . 
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Proof. We have 

I 9R a [l}f + Ig R [2]g = [l][f(b, c) - f(b, d)} + e IS[2][g(b, d) - g(a, d)} . 
By equation (@}, 

f(b, c) - f(b, d) = -- J2 IW [2}f(b, d) , 
d) - g(a, d) = - ; a I?[l] a D%[l]g(b, d) . 

Therefore, 

ImW + ISr[2]9 = -jm^cI2mcD2l2]f(b,d) + C I%[2}±- I?[l] a DZ[l]g(b,d) 

□ 

Theorem 2 (Fractional Gauss's theorem for a parallelepiped). Let F — (F x ,F y ,F z ) be a contin- 
uous vector field in a domain that contains W = [a,b] x [c,d] x [e, /]. If the boundary of W is a 
closed surface dW , then 

{I$ w ,F) = -I^Div^F. (5) 

Proof. The result follows by direct transformations: 

(Idw F ) = ^w[ 2 ' 3 ] f i +^H'[ 1 ' 3 l- f, j +^9w[l>2]-F« 

= c I^}eIfm(F x (b,dJ) -F x (a,d,f)) + a I?[l] e I?[3](F y (b,d,f) -F y (b,c,f)) 

+ a I? [l]J2 [2] (F z (b, d, f) - F z (b, d, e)) 

= ^ a ini}ci2meifm( a Dnm(b, d, /) + c m^F y (b, n + eD^m^, d, 

= ±I&( a DZ[l]F x + c D^[2]F y + e DJ[3]F z ) 

□ 

Let S be an open, oriented, and nonintersecting surface, bounded by a simple and closed curve 
dS. Let F = [F x ,F y ,F z ] be a continuous vector field. Divide up S by sectionally curves into N 
subregions Si, Sa, • ■ ., Sn- Assume that for small enough subregions each Sj can be approximated 
by a plane rectangle Aj bounded by curves C\, C2, ■ ■ ., Cn- Apply Green's theorem to each 
individual rectangle Aj. Then, summing over the subregions, 

3 3 

Furthermore, letting N — > 00 

Y,^(V%><F)^±(I%,Cvrl<* s F) 

3 

while 

3 

We conclude with the fractional Stokes formula 

l(^,Curl^) = ^ 5 F. 
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4 Fractional calculus of variations with multiple integrals 

Consider a function w — w(x, y) with two variables. Assume that the domain of w contains the 
rectangle R = [a, b] x [c, d] and that w is continuous on R. We introduce the variational functional 
defined by 



J(w) = I%L (x, y, w(x, y), a D a x [l]w(x, y), c D«[2]w(x, y)) 

b r d (6) 
- 1 1 L(x,y,w, a D^[l}w, c D^[2}w)(b-xr- 1 (d^y) a - 1 dydx 



:= a 

We assume that the lagrangian L is at least of class C 1 . Observe that, using the notation of the 
(dt) a integral as presented in [28], (|6]) can be written as 

JH = f [ d L{x,y,w(x,y), a D«[l}w(x,y), c D«{2}w(x,y)) {dy) a {dx) a . (7) 

J a J c 

Consider the following FCV problem, which we address as problem (P). 

Problem (P): minimize (or maximize) functional J defined by ([7]) with respect to the set of 
continuous functions w(x,y) such that w\qr = ip(x,y) for some given function tp. 

The continuous functions w(x, y) that assume the prescribed values w\qr = <p{x, y) at all points 
of the boundary curve of R are said to be admissible. In order to prove necessary optimality 
conditions for problem (P) we use a two dimensional analogue of fractional integration by parts. 
Lemma [1] provides the necessary fractional rule. 

Lemma 1. Let F , G, and h be continuous functions whose domains contain R. If h = on dR, 
then 

r [G(x, y) a D a x [l]h{x 1 y) - F(x, y) c D"[2]h(x, y)}(b - x) a -\d - y)^ 1 dydx 



u 

i, ,..,/ 

v 



[( D«[l]G(i, y) - cD«{2]F(x, y))h(x, y)](b - x) a -\d - y)"" 1 dydx . 



Proof. By choosing / = F ■ h and g = G ■ h in Green's formula, we obtain 

ImlM(Fh) + I% R [2](Gh) = ijg[ D?[l]G ■ h + G ■ a D^[l]h - C D%[2]F -h-F- c D<* d [2]h] , 
which is equivalent to 

i/£ [G ■ a D%[l]h - F ■ c D a d [2]h] = I$ R [l](Fh) + I§ R [2](Gh) - [{ a D%[l]G - c D a d [2]F)h] . 
In addition, since h = on dR, we deduce that 

iS [G ■ aD£[l]h - F ■ c D a d [2]h] = -1% [{ a D%[l]G - c D a d [2]F)h] . 
The lemma is proved. □ 

Theorem 3 (Fractional Euler-Lagrange equation). Let w be a solution to problem (P). Then w 
is a solution of the fractional partial differential equation 

d 3 L- a D«[l]d4L- c D%[2]d 5 L = Q, (8) 

where by diL, i = 1, . . . , 5, we denote the usual partial derivative of L(-, •, •, •, •) with respect to its 
i-th argument. 
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Proof. Let He a continuous function on R such that ft e on dR, and consider an admissible 
variation w + eh, for e taking values on a sufficient small neighborhood of zero. Let 

j(e) = J{w + eh). 

Then j'(0) = 0, i.e., 

pb pd 

a 2 (d 3 Lh + d 4 L a D«[l}h + d 5 L c D*[2}h){b-x) a - 1 (d~y) a - 1 dydx = Q. 

J a J c 

Using Lemma [T] we obtain 

pb pd 

a 2 (d 3 L- a D«ll}d 4 L- c D«l2}d 5 L)h(b-x) a - 1 (d~y) a - 1 dydx = 0. 

J a J c 

Since ft is an arbitrary function, by the fundamental lemma of calculus of variations we deduce 
equation ©. □ 

Let us consider now the situation where we do not impose admissible functions w to be of fixed 
values on OR. 

Problem (P'): minimize (or maximize) J among the set of all continuous curves w whose 
domain contains R. 

Theorem 4 (Fractional natural boundary conditions). Let w be a solution to problem (P'). Then 
w is a solution of the fractional differential equation (jHJ) and satisfies the following equations: 

1. d 4 L(a,y,w(a,y), a D*[l}w(a,y), c D%[2}w(a,y)) = for all y G [c,d]; 

2. 8 4 L(b, y, w(b, y), a D%[l]w(b, y), C D% [2]w(b, y)) = for all y e [c, d}; 

3. d 5 L(x,c,w(x,c), a D%[i\w(x,c), c D*[2]w(x,c)) = for all x e [a,b]; 

4. d 5 L(x,d,w(x,d) ta D%[l]w(z,d), c D%[2]w(x,d)) = for all x 6 [a, b]. 
Proof. Proceeding as in the proof of Theorem [3] (see also Lemma [T]), we obtain 

= a 2 / (d 3 Lh + d 4 L a D«[l]h + d 5 L c D*[2}h) (b - x)*- 1 {d - y)^ 1 dydx 

J a J c 

pb pd 

= a 2 / (d 3 L- a D^[l}d 4 L- c D^[2}d 5 L)h(b~x) a - 1 (d-y) a ~ 1 dydx 

J a J c 

+ a\I§ R [2](d 4 Lh) - a\I§ R [l](d 5 Lh) , (9) 

where h is an arbitrary continuous function. In particular the above equation holds for ft e on 
dR. For such ft the second member of © vanishes and by the fundamental lemma of the calculus 
of variations we deduce equation ((SJ). With this result, equation © takes the form 

= ^ d 4 L(b, y, w(b, y), a DZ[l]w(b, y), c D«[2]w(b, y)) h(b, y)(d - y^dy 

d 



d 4 L(a, y, w(a, y), a D%[l]w(a, y), C D* [2}w(a, y)) h(a, y){d - y) a 1 dy 

b 

- I d 5 L(x,c,w(x,c), a D%[l]w(x,c), c D*[2}w(x,c)) h(x,c)(b - x) a ^dx 



(10) 



d 5 L(x, d, w(x, d), a D*[l}w{x, d), c D°[ [2)w(x, d)) h(x, d)(b - x^dx 
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Since h is an arbitrary function, we can consider the subclass of functions for which h = on 

[a, b] x {c} U [a, b] x {rf} U {6} x [c, d]. 
For such h equation (fit))) reduce to 

= y a 4 i(o, y, w(a, y), a D%[l]w(a, y), c D%[2]w(a, yj) h(a, y)(d - yT^dy. 

By the fundamental lemma of calculus of variations, we obtain 

d i L(a,y,w(a,y), a D%[l]w(a,y), c D™[2}w(a,y)) = for all y e [c,d] . 

The other natural boundary conditions are proved similarly, by appropriate choices of h. □ 

We can generalize Lemma Q] and Theorem |3] to the three dimensional case in the following way. 

Lemma 2. Let A, B, C , and rj be continuous functions whose domains contain the parallelepiped 
W. Ifr) = on dW, then 

Ifr(A ■ a D%[l]r, + B ■ C D%[2} V + C ■ e D?[3} V ) = -I%{[ a D%[l]A + e D%[2}B + e DJ[3}C] r,) . (11) 

Proof. By choosing F x = r]A, F y = nB, and F z = rjC in ([5]), we obtain the three dimensional 
analogue of integrating by parts: 

(A ■ a Dt [1] V + B ■ c Bl [2}r, + C ■ e DJ [3}r,) 

= -I^{[ a D%[l]A + C D%[2]B + e Df[3]C] r,) + a\(Ig w , [r>A, V B, V C}). 

In addition, if we assume that rj = on dW, we have formula (fTTj) . □ 

Theorem 5 (Fractional Euler-Lagrange equation for triple integrals). Let w — w(x,y,z) be a 
continuous function whose domain contains W — [a,b] x [c,d] x [e, /]. Consider the functional 

J(w) = IfrL (x, y, z, w(x, y, z), a D%[l]w(x, y, z), c L>^[2]w(a;, y, z), e D^[3]w{x, y, z)) 
= {( [ L( X ,y,z,w, a D%[l]w, c D%[2}w, e D%{3]w) (dz) a (dy) a (dx) a 

J a J c J e 

defined on the set of continuous curves such that their values on dW take prescribed values. Let 
L be at least of class C 1 . If w is a minimizer ( or maximizer) of J, then w satisfies the fractional 
partial differential equation 

d 4 L - a DZ{l]d 5 L - c D%[2]d e L - e DJ[3]d 7 L = 0. 

Proof. A proof can be done similarly to the proof of Theorem [3] where instead of using Lemma [1] 
we apply Lemma [5J □ 

5 Applications and possible extensions 

In classical mechanics, functionals that depend on functions of two or more variables arise in a 
natural way, e.g., in mechanical problems involving systems with infinitely many degrees of freedom 
(string, membranes, etc.). Let us consider a flexible elastic string stretched under constant tension 
t along the x axis with its end points fixed at x — and x = L. Let us denote the transverse 
displacement of the particle at time t, t\ < t < whose equilibrium position is characterized 
by its distance x from the end of the string at x — by the function w = w(x,t). Thus w(x,t), 
with < x < L, describes the shape of the string during the course of the vibration. Assume a 
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distribution of mass along the string of density a — a(x). Then the function that describes the 
actual motion of the string is one which renders 



an extremum with respect to functions w(x,t) which describe the actual configuration at t = t\ 
and t = <2 and which vanish, for all t, at x = and x — L (see [43J, p. 95] for more details). 

We discuss the description of the motion of the string within the framework of the fractional 
differential calculus. One may assume that, due to some constraints of physical nature, the 
dynamics do not depend on the usual partial derivatives but on some fractional derivatives o-D" [l]w 
and tl D"[2]w. For example, we can assume that there is some coarse graining phenomenon - see 
details in 26,29. In this condition, one is entitled to assume again that the actual motion of the 
system, according to the principle of Hamilton, is such as to render the action function 



where R = [0, L] x [t\, £2], an extremum. Note that we recover the classical problem of the vibrating 
string when a — > 1~. Applying Theorem [3] we obtain the fractional equation of motion for the 
vibrating string: 



This equation becomes the classical equation of the vibrating string (cf., e.g., j43j p. 97]) if a 1 




We remark that the fractional operators are non-local, therefore they are suitable for con- 
structing models possessing memory effect. In the above example, we discussed the application of 
the fractional differential calculus to the vibrating string. We started with a variational formula- 
tion of the physical process in which we modify the Lagrangian density by replacing integer order 
derivatives with fractional ones. Then the action integral in the sense of Hamilton was minimized 
and the governing equation of the physical process was obtained in terms of fractional derivatives. 
Similarly, many others physical fields can be derived from a suitably defined action functional. 
This gives several possible applications of the fractional calculus of variations with multiple inte- 
grals as was introduced in this paper, e.g., in describing non-local properties of physical systems 
in mechanics (see, e.g., [TTJ [T^l [321 13S1 SI]) or electrodynamics (see, e.g., [TU1 132] ). 

We end with some open problems for further investigations. It has been recognized that 
fractional calculus is useful in the study of scaling in physical systems [131115]. In particular, there 
is a direct connection between local fractional differentiability properties and the dimensions of 
Holder exponents of nowhere differentiable functions, which provide a powerful tool to analyze 
the behavior of irregular signals and functions on a fractal set [H [33J . Fractional calculus appear 
naturally, e.g., when working with fractal sets and coarse-graining spaces [29[ [30], and fractal 
patterns of deformation and vibration in porous media and heterogeneous materials [2]. The 
importance of vibrating strings to the fractional calculus has been given in [18] , where it is shown 
that a fractional Brownian motion can be identified with a string. The usefulness of our fractional 
theory of the calculus of variations with multiple integrals in physics, to deal with fractal and 
coarse-graining spaces, porous media, and Brownian motions, are questions to be studied. It 
should be possible to prove the theorems obtained in this work for a general form of domains and 
boundaries; and to develop a fractional calculus of variations with multiple integrals in terms of 
other type of fractional operators. An interesting open question consists to generalize the fractional 
Noether-type theorems obtained in ETJ [22] to the case of several independent variables. 
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